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The aim of this study was to analyse the relationships between polymorphisms in four candidate genes (AR, 
DATl, DRD2, and COMT) and aggression in men from a traditional society of East African pastoralists, the 
Datoga. Buss and Perry's Aggression Questionnaire was used to measure aggression. The number of CAG 
repeats in the AR gene was negatively correlated with physical aggression, anger, and hostility. Among the 
genes of the dopaminergic system, a significant single-gene effect was detected only for DRD2 with regard to 
anger. At the level of a two-gene model, a significant effect for DRD2 and a tendency for DATl were 
observed for the DAT1-DRD2 gene pair regarding hostility, and two tendencies were observed for the 
interaction effect of the DA Tl -COMTpair regarding anger and hostility. These data suggest a probable link 
between physical aggression and direct fitness caused by strong sexual selection in Datoga men. 



Aggressive behaviour in humans is highly sexually dimorphic, with men exhibiting more violent (physical) 
same-sex aggression than women in almost all cultures 14 , which is most likely attributable to the greater 
impulsiveness of men and a stronger fear of physical danger in women 5,6 . A meta-analysis of 24 genetic 
studies on aggression demonstrated that heritability accounts for approximately 50% of the variance in aggres- 
sion 7,8 , with the heritability of aggressive behaviour being higher in men than in women 9 . Aggressive behaviour 
(physical and verbal) and its negative emotional bases (anger and hostility) are the products of complex inter- 
actions between many genes associated with the secretion of steroid hormones and the sensitivity of tissues to 
these hormones. These genes include those associated with the dopaminergic and serotonergic systems (both 
receptors and transmitters) and likely many others 1012 . 

The results of some studies suggest a relationship between exposure to testosterone in utero and subsequent 
aggression in adulthood as a result of masculinisation of the brain 13,14 . Individual behavioural reactions to 
testosterone have also shown a close dependence on the sensitivity of brain tissue to this steroid. Androgen 
receptor genes are located on the X chromosome and are therefore only inherited by men from their mothers. 
Earlier studies, predominantly conducted among representatives of industrial populations, appear to support the 
hypothesis of a negative correlation between an increasing the number of CAG repeats in the androgen receptor 
gene (AR) and physical aggression and/or the 2D : 4D ratio (the ratio between the lengths of the second and fourth 
digits, an indicator of the masculinisation of the foetus in utero) 15 . Another study reported that fewer CAG repeats 
in AR were more frequently found in male rapists and murderers than in control subjects 16 . Indeed, an increased 
number of these repeats in AR results in receptors with lower androgen sensitivity 17,18 . Manning and colleagues 
demonstrated in a male sample from an industrial population that masculine finger ratios are associated with 
androgen receptor alleles with fewer CAG microsatellite repeats in the terminal domain 15 . In contrast, a recent 
study performed by our group in traditional hunter-gatherers, the Hadza of Tanzania, revealed no correlation 
between the number of CAG repeats and aggression or the right-hand 2D : 4D ratio in Hadza men 1 '. Thus, the 
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recent data on the associations between AR polymorphisms and 
aggression are somewhat contradictory and demand a more thor- 
ough investigation in other small-scale societies. 

Many recent molecular genetic studies on aggression have focused 
on the genes involved in dopamine neurotransmission. The distribu- 
tions of polymorphic loci in the dopamine receptor genes DRD4 and 
DRD2 and the dopamine transporter gene DAT1 were analysed to 
identify an association between specific alleles and aggressive beha- 
viour (particularly obsessive-compulsive disorder, hostility, and 
physical aggression), some mental disorders, and alcohol addic- 
tion 20-27 . Another gene associated with aggression is the gene encod- 
ing catechol-O-methyltransferase (COMT), which plays a key role in 
dopamine catabolism in the prefrontal cortex: the V158M poly- 
morphism in the fourth exon of the COMT gene leads to a 40% 
reduction of enzymatic activity in carriers of the M allele 28 . 

Dopamine plays a critical role in reward processing and is a potent 
neuromodulator of ventral striatum (VS) reactivity, which is known 
to be highly involved in reward processing. VS reward reactivity may 
be a key neurobiological pathway through which DAT1 and DRD2 
polymorphisms contribute to the variability of behavioural impul- 
sivity 29 . 

DAT1 and DRD2 polymorphisms were examined in adolescents 
with pathological aggression, and it was found that carriers of the 
DAT1 10/10 genotype were overrepresented in a sample of aggres- 
sors compared to the control group and that DRD2 Al/Al and All 
A2 were also overrepresented in aggressors 30 . Vaughn and collea- 
gues 31 suggested DAT1 as one of the main genetic markers of anti- 
social behaviour in adolescent criminals. According to another study, 
the DAT1 10/10 and 10/9 genotypes are present in violently delin- 
quent young adults twice as often as in controls 26 ' 32 " 34 . Guo and co- 
authors 26 demonstrated that heterozygotes for the DRD2 A1/A2 
variant presented higher aggression scores than Al/Al or A2/A2 
homozygotes; these authors further suggested that there is no cor- 
relation or epistasis between DRD2 and DAT1. 

Remarkably, the data from a nationally representative sample of 
American youth revealed that genetic risk factors measured based on 
two dopamine polymorphisms (DAT1 and DRD2) confound the 
association between paternal incarceration and a child's arrest 35 . 

A serious problem with the current studies on the molecular gen- 
etics of aggression and antisocial behaviour is that such research 
frequently cannot be replicated 36 . Moreover, a recently conducted 
meta-analysis of genetic association studies on violence and aggres- 
sion provided a rather pessimistic outlook, indicating that the can- 
didate gene approach has not succeeded in identifying genes 
associated with these outcomes 37 . One of the general explanations 
for this difficulty is that "there are few, if any, loci with large effect 
size, and it is becoming increasingly obvious that it will be necessary 
to consider the impact of genes, not in isolation, but as part of a 
multifactorial miasma including both other genetic factors and the 
environment" 9 . Because most studies are based on clinical samples or 
samples that are not nationally representative and each sample is 
likely differentially exposed to particular environmental conditions, 
certain genes may only have effects when they are paired with certain 
environments 38 . Another possible reason for the lack of reproducib- 
ility in this arena is population stratification. Indeed, the allelic fre- 
quencies of some genes vary significantly across different racial/ 
ethnic categories 39 , and genetic effects may vary across racial/ethnic 
groups because of (a) differences in allelic distributions and (b) dif- 
ferences in exposure to environmental risk factors 38 . 

Thus, there is an urgent need to examine the associations between 
certain genes and aggressive behaviour in samples of various racial/ 
ethnic groups. In this study, we examined the single-gene effects of 
the AR gene and three dopaminergic genes (DAT1, DRD2, and 
COMT) as well as the two-gene associations of components of the 
dopaminergic system with aggression in a sample of Datoga, who are 
traditional East African pastoralists. 
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The following hypotheses were tested: 

1. Men showing lower AR CAG repeat numbers exhibit higher 
aggression scores, particularly for physical aggression and 
negative emotions. 

2. Male carriers of the +10 type of DAT1 (homozygous or het- 
erozygous) present higher aggression scores. 

3. Men heterozygous for DRD2 display higher aggression scores 
compared to both types of homozygotes. 

4. Men homozygous for the high-activity Val/Val genotype of 
COMT show higher aggression scores. 

5. Two-gene interaction effects between the three genes in the 
dopaminergic system (DAT1, DRD2, and COMT) provide 
important information on the associations of these genes with 
various types of aggression. 

Results 

The means and standard deviations for Buss and Perry's Aggression 
Questionnaire subscales were computed as well as minimum, and 
maximum scores for Datoga males (Table 1). 

We performed a simple regression analysis to test the single-gene 
effect of the AR gene on aggression and detected significant effects of 
the AR CAG polymorphism on physical aggression and hostility, 
with carriers of smaller numbers of CAG repeats showing a higher 
rating (Tables 2, and Figs. la,b,c). The results for AR and the self- 
rated scores for anger exhibited a trend in the same direction 
(Table 2). Therefore, our data show that men with fewer CAG repeats 
behave more aggressively than their fellow males with a greater 
number of AR CAG repeats. 

We split our sample according to DAT1 genotype, as suggested 
previously by a number of authors 26,30,32 . The first group included 
men with the 9/9 genotype of the DAT1 ( — 10) type, and the second 
group included individuals with the 9/10 or 10/10 genotype of the 
DAT1 (+10) type. One-way ANOVA with DAT1 type as the inde- 
pendent variable and aggression traits as the dependent variables 
revealed no significant association between the +10 type and all 
AQ scales (Table 3). 



Table 2 | Simple regression analyses. Main effects of ARgene poly- 
morphism on physical aggression, anger and hostility in Datoga 
males 

AQ Scales N R 2 (3 t P 

PhA 87 0.105 -0.3234 -3.0189 0.0042** 

Ang 85 0.080 -0.2827 -2.6025 0.01 97 + 

Hostil 86 0.069 -0.2531 -2.3104 0.0142* 



PhA - physical aggression, Ang - anger, Hostil - hostility, N - the sample size, R 2 - proportion of 
the total variability explained by the factor effect, p - coefficient of regression, t- f-statistic; 
*P<0.05; 
**P<0.01; 

+ P < 0.1 (adjusted for multiple correlated outcomes). 
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Figure 1 | The distributions of the physical aggression (A), anger (B), and hostility (C) subscales as correlated with the number of CAG repeats in the 
AR gene. The ordinate numbers designate the self-rated scores according to the subscales of Buss and Perry's Aggression Questionnaire. 
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Our sample was also split into three groups based on DRD2 poly- 
morphisms, with group 1 including men with the Al/Al genotype, 
group 2 consisting of men with the A1/A2 genotype, and group 3 
comprising individuals with the A2/A2 genotype, following the divi- 
sions originally employed in a previous study by Guo and collea- 
gues 26 . Significant single-gene effects of DRD2 in the Datoga men 
were found only with regard to anger (Fig. 2). Lastly, we divided our 
sample into three groups according to COMT polymorphisms: group 
1 (V/V), group 2 (V/M), and group 3 (M/M); however, no single- 
gene effect of COMT was identified in this study (Table 3). 

We subsequently performed a factorial two-way ANOVA to test 
both the main and interaction effects of the three genes (J) ATI, 
DRD2, and COMT) on three aggression subscales (Tables 4). The 
two-gene model revealed two significant main effects of the DRD2 
gene on anger and hostility, showing a tendency for a main effect of 
DAT1 on hostility. Tendencies of DA Tl X COMT interaction effects 
on anger and hostility, but not on physical aggression, were also 
observed (Fig. 3a,b). 

Discussion 

Overall, we found that the polymorphisms in all four genes tested in 
this study affected the self-rated aggression scores of the Datoga men, 
with either single-gene effects or interaction effects between these 
genes being detected. Our data also confirmed our initial hypothesis 
regarding the single-gene effect of AR on the self- rated scores for the 
three AQ scales. Our results regarding the association between 
aggression and the number of CAG repeats in AR were in the pre- 
dicted direction: the levels of physical aggression, anger, and hostility 
among the Datoga men were higher in individuals with fewer CAG 
repeats. Interestingly, we recently demonstrated that individuals 
with a small number of CAG repeats are overrepresented among 
world-class judo sportsmen compared to the general Russian popu- 
lation 40 . At the same time, these data contradict those obtained by 
Hurd, Vaillancourt, and Dinsdale 41 in Canadian students, where men 
with more sensitive androgen receptors tended to score lower for 
physical aggression and anger. The results also contradict our early 
findings in Hadza men, for whom no correlations were found 
between the AQ subscales and the number of AR CAG repeats 42 . A 
reasonable explanation for these differences may be the existence of 
gene-environment effects and may be particularly attributable to 
differences in cultural attitudes toward aggression among the Hadza 
and Datoga populations. 



The significant single-gene effects of DAT1 polymorphisms on 
anger and hostility were in the direction predicted in the second 
hypothesis. Therefore, not only is +10 DAT1 overrepresented in 
samples of antisocial adolescents in the USA 31 , but it is possible that 
these genotypes are positively selected in pastoralist cultures because 
such personalities are better suited to both cattle raiding and cattle 
protection. Although significant single-gene effects of the DRD2 
polymorphism on anger and hostility were also demonstrated, 
contrary to our expectations, the carriers of the heterozygous type 
displayed the same high levels of these characteristics as the homo- 
zygous A2/A2 men. We did not find any single-gene effect of COMT 
in this study; therefore, our fourth hypothesis should be rejected. 
Lastly, we demonstrated tendencies toward interaction effects 
between DAT1 and COMT with regard to anger and hostility in 
Datoga males. According to our findings, these interaction effects 
are only essential for carriers of the different COMT genotypes and 
the -10 DAT1 type. 

Guo and colleagues 26 demonstrated that an interaction between 
DAT1 and DRD2 produces essentially the same effect on aggression 



< 




A1/A1 A1/A2&A2/A2 
DRD2 Type 

Figure 2 | The distribution of the anger subscale as correlated with the 
DRD2type (1 - Al/Al and2 - A1/A2 and A2/A2). The designations of the 
numbers are described in Figure 1 . 
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observed in single-gene models, which was confirmed in our analysis 
of anger and hostility. 

However, we did not test putative associations of AR gene poly- 
morphisms with the three genes involved in the dopaminergic sys- 
tem, as there is still no evidence of their functional relationship. 

Because these data were collected in a single population and due to 
the highly similar cultural and social environments of the respon- 
dents, we believe that we were able to minimise the inter- individual 
variation in the environmental influence. Moreover, only males were 
included in our analysis, as there is convincing evidence to support 
significant sex-based differences in gene-environment interac- 
tions 43 . The Datoga are polygynous pastoralists, and the trajectories 
for achieving higher reproductive success are radically different in 
men and women. In contrast to the egalitarian and relatively peaceful 
Hadza hunter-gatherers, the men of the Datoga, another indigenous 
culture in the same region 44,45 , have been selected for their aggression, 
particularly their physical aggression, as have other African pastor- 
alists 46,47 . Indeed, men with fierce characters and warrior skills may 
be more successful in their mating and parental efforts because they 
can acquire and protect cattle, consequently showing higher repro- 
ductive success (more wives and children) 48 . Furthermore, our data 
suggest a possible link between physical aggression and direct fitness 
caused by strong sexual selection in Datoga men, while such a rela- 
tionship most likely no longer exists in Westernised societies. 

Of course, the present work has a number of limitations. First, the 
sample size was not large. However, because all our respondents were 



from the same tribal community, practiced traditional rites, partici- 
pated in the traditional economy, and lived in the same area, the 
inter- individual variability in this respect was somewhat limited. 
As mentioned above, previous findings regarding the association 
between genes and particular behavioural traits are highly contra- 
dictory 36-37 . One possible reason for these disparities maybe rooted in 
environmental factors 38 , and the present data together with our data 
on the Hadza 42 illustrate a situation favouring this suggestion. 
Additionally, particular polymorphic loci may account for a rela- 
tively small proportion of phenotypic variance, and these effects 
may not be significant in relatively small samples of individuals with 
different socialisation experiences. If more than two genes/poly- 
morphisms and many behavioural traits are selected for analysis, 
the necessity of making corrections for multiple comparisons is 
problematic. One possible solution for this situation is to obtain a 
very large sample; another possibility is the generation of biologically 
founded multilocus genetic profiles reflecting the cumulative effects 
of multiple polymorphic loci of known functionality on a specific 
behaviour 49 . In accord with this approach, Beaver and Chavian cre- 
ated the Genetic Risk Index for criminal involvement in Hispanics 
based on three dopaminergic gene polymorphisms and demon- 
strated its effectiveness 38 . Other authors proposed the construction 
of a multilocus genetic profile score based on five functional poly- 
morphic loci for dopamine signalling and reported that it accounted 
for 10.9% of the inter- individual variability in reward-related VS 
reactivity 50 . The dopamine (DA) index was found to be associated 
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Figure 3 | The distribution of the anger (A) and hostility (B) as correlated 
with the DAT1 x COMTeffect. The designations of the ordinate numbers 
are described in Figure 1. 

with ventral striatum activity in an experimental study and was 
recently applied in a psychiatric study. However, we decided not to 
use this index in our work because we do not know how it may be 
associated with the activity of other brain structures that are activated 
during aggression and their control in humans, particularly the pre- 
frontal zone 51 . We refrained from applying such an approach in this 
study because (i) the initial assumption that particular polymorph- 
isms act additively and not interactively to influence dopamine sig- 
nalling should be tested in more detail, and (ii) more evidence is 
needed that all of the loci used to generate the multilocus genetic 
profile score are truly invested equally in behaviour. 



Methods 

Ethics statement. Institutional approvals, including those from the University 
(Moscow State University Ethics Committee) and local governmental agencies 
(including the Tanzanian Commission for Science and Technology), were obtained 
prior to conducting this study. All subjects gave their informed, verbal consent prior 
to participation. Verbal consent was deemed appropriate given the low literacy rates 
among traditional Datoga and was specifically approved by the University EC and 
Tanzanian agency. 

Study sample and data collection. The data were collected during field studies from 
2006 to 2012 in the United Republic of Tanzania. The sample analysed here includes 
138 adult Datoga men, with a mean age of 34.09 ± 12.12 years (median of 33.5 years; 
range, 18-70 years). 

The participants were personally interviewed to determine their age, sex, ethnicity, 
and personal history, as described previously 42 . All participants gave their consent 
prior to participating, and the average pay per subject was equal to US$2. All inter- 
views were conducted in a one-to-one dialogue with the respondents. All questions 
were read aloud and explained if necessary. Many of the Datoga men commented in 
these interviews that aggression is an appropriate way of 'saving face', for both 
themselves and their relatives, and of protecting their property. When asked about 
family violence, all the men mentioned that they were physically punished during 
childhood by their parents (father, mother, or both). Self-rated aggression scores were 
obtained using Buss and Perry's Aggression Questionnaire (AQ) 52 , which consists of 

29 statements, grouped into four subscales: physical aggression (nine items), verbal 
aggression (five items), anger (seven items), and hostility (eight items). The first two 
subscales, physical and verbal aggression, represent direct aggression, and the other 
two subscales pertain to the individual's emotional basis for aggressive behaviour. AQ 
uses a Likert scale, ranging from 1 (extremely uncharacteristic) to 5 (extremely 
characteristic). This tool was employed because the original English version of the AQ 
has shown moderate construct validity and high test-retest reliability 53 and has been 
widely applied in aggression studies. In the present work, we used the ki-Swahili 
version of the AQ, which was initially administered by us and was previously used in 
our research on the Hadza 42 . For the details of the data collection procedure, see 
Butovskaya et al 42 . Here, we present data for three subscales of the AQ: physical 
aggression, anger, and hostility. Cronbach's Alpha values for these behavioural 
measures were as follows: physical aggression, 0.64; anger, 0.63; and hostility, 0.66. 

Although the main conclusions drawn in this work regarding the subjects' 
aggression profiles were based on self-reported scores, rather than on direct obser- 
vations, we do not think that this method of data collection significantly distorted the 
results. A meta-analysis performed earlier by Book and colleagues 54 demonstrated 
that "the measure of aggression (behavioural or self-report) did not have any effect on 
the relationship between testosterone and aggression," and there is no reason to 
assume that this situation differs with regard to genetic data. 

Genotyping. Buccal epithelium samples were collected for DNA analysis during the 
field surveys, and DNA was later isolated using the DIAtom 1M DNA Prep 200 
extraction kit (IsoGene Lab, Moscow, Russia), which is intended for the isolation of 
DNA from various biological materials, according to the manufacturer's protocol. 

AR CAG repeats, DRD2 TaqlA, the DA Tl VNTR 3 ' -noncoding region and COMT 
Val-158-Met loci polymorphisms were typed via polymerase chain reaction (PCR), 
which was performed with the GenePak® PCR MasterMix Core (IsoGene Lab, 
Moscow, Russia) reagent kit, according to the recommendations of the manufac- 
turer 55 . 

Amplification of the target loci was performed using an MJ Research PTC- 100 
thermocycler with the following program: 1 cycle of 1 min at 94°C for denaturation, 

30 sec at X°C for annealing, and 30 sec at 72 IJ C for synthesis. This was followed by 30 
cycles of 30 sec at 94 C, 30 sec at X C, and 15 sec at 72°C, with a final elongation at 
72 'C for 3 min. X denotes the annealing temperatures of individual primer pairs 
(Table 5). 

The reaction products for the AR locus were analysed using an ABI PRISM 3100- 
Avant automated DNA sequencer (Applied Biosystems, Foster City, USA). 

To identify the TaqI polymorphism (T/C), the DRD2 amplification product was 
subdivided into 10 ul aliquots, one of which was treated with the TaqI restriction 
endonuclease overnight at 65 C. For the COMT locus, the PCR products were also 



Table 5 | Summary of 


primers for the PCR amplification of the four loci studied 




Locus 


Primer Sequences 


Annealing Temperatures 


AR(CAG)„ 


F* : 5'-(FAM)-TCCAGAGCGTGCGCGAAGTGAT-3 ' 


56-C 




R: 5 '-(FAM)-CGACTGCGGCTGTGAAGGTTG-3 ' 




DRD2-Taq 1 A 


F: 5 '-CCGTCGACGGCTGGCCAAGTTGTCTA-3 ' 


68°C 


R : 5 'GCGTCGACCCTTCCTGAGTGTC ATC A3 ' 




DAT1 VNTR 


F: 5 '-TGCGGTGTAGGGAACGGCCTGAG-3 ' 


68°C 




R: 5 '-CTTCCTGGAGGTCACGGCTCAAGG-3 ' 




COMTM ,58 V 


F: 5 '-TACTGTGGCTACTCAGCTGTGC-3 ' 


55°C 




R: 5 'GTGAACTGTGTGTGAACACC-3 ' 




*F - forward primer, R - reverse primer. 
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subdivided into 10 ul aliquots, and one of these aliquots was treated with Nlalll 
overnight at 37 °C, 

The amplification and restriction products were separated via electrophoresis in a 
2% (DRD2 and DAT1) or 3% (COMT) agarose gel and stained with ethidium 
bromide. The results were photographed and analysed with a BioDocAnalyze device 
(Biometra, Goettingen, Germany). 

Statistics. Statistical treatment of the data was performed with SPSS 13.0 for 
Windows (IBM, Armonk, NY, USA). The means and standard deviations were 
computed for the studied traits (Table 1). The Kolmogorov-Smirnov test was 
employed to test whether our behavioural data were normally distributed. As all of the 
obtained probability values were greater than 0.1, parametric methods of statistical 
analysis were suitable for evaluating the dependence of the aggression subscales on 
the genotypes. The scores at all four aggression subscales were positively correlated 
with one another (with r values ranging from 0.386 to 0.497); thus, we used the 
correction for the critical p value for multiple correlated outcomes {http://gump.qimr. 
edu.au/general/daleN/matSpD/) previously applied by Hurd, Vaillancourt, and 
Dinsdale 41 . The adjusted critical p value for all of the aggression subscales was 
0.01695. Because we tested the association of the AR gene and a set of three 
dopaminergic genes on aggression scales separately, the following corrections of the 
critical p-values were suggested. The adjusted critical value for AR was set to p < 
0.017/l(gene) — 0.017 (two-tailed), whereas p ^ 0.034 was regarded as a trend. The 
significance levels for the three dopaminergic genes were calculated as p < 0.017/ 
3(genes) — 0.0057 (two-tailed); a level of probability ranging from 0.0057 to 0.0113 
was considered a tendency. 

A linear regression was applied to analyse the associations of the three aggression 
subscales with the number of CAG repeats in the AR gene. One-way and two-way 
ANOVAs were conducted to determine the effects of the three candidate dopamine 
genes on the variations in the examined behavioural traits. We did not test the 
interaction effects of AR and the three dopamine genes in association with aggression 
because they represent different functional systems. 

A genetic population analysis revealed agreement with Hardy- Weinberg equilib- 
rium for the genotypes of the DAT1 and COMT genes (Hobs = 0.505/Hexp = 0.505, 
p = 0.835, and Hobs = 0.491/Hexp = 0.443, p = 0.298, respectively), but not for the 
DRD2 gene (Hobs = 0.692/Hexp = 0.495, p <C 0.001, with a corrected level of 
significance of a/3 - 0.05/3 - 0.017). 
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